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We report a search for diphoton events with large missing transverse energy produced in pp collisions at 
•Js = 1 .96 TeV. The data were collected with the DO detector at the Fermilab Tevatron Collider, and correspond 
to 6.3 fb _1 of integrated luminosity. The observed missing transverse energy distribution is well described 
by the standard model prediction, and 95% C.L. limits are derived on two realizations of theories beyond the 
standard model. In a gauge mediated supersymmetry breaking scenario, the breaking scale A is excluded for 
A < 124 TeV. In a universal extra dimension model including gravitational decays, the compactification radius 
R c is excluded for R' 1 < All GeV. 
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In the standard model (SM), events with two high trans- 
verse momentum photons (yy) and large missing transverse 
energy (fir) are produced at a small rate in pp collisions. 
This final state is therefore sensitive to contributions from pro- 
cesses beyond the SM (BSM). We report a search for yy events 
with large Tf.j produced in pp collisions recorded using the DO 
detector at the Fermilab Tevatron Collider. The sensitivity is 
assessed for two benchmark BSM models, gauge mediated 
supersymmetry (SUSY) breaking (GMSB) (T) and universal 
extra dimensions (UED) 

In GMSB models, the masses of the SUSY partners to SM 
particles arise from SM gauge interactions and are propor- 
tional to the effective SUSY breaking scale A. As the grav- 
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itino (G) does not participate in SM gauge interactions, it has 
a small mass [3 1 and is the lightest SUSY particle. Assuming 
R parity conservation [4|, the SUSY process with the largest 
cross section at the Tevatron would be chargino and neutralino 
pair production (%f %?> X\X\) 0' followed by decay chains 
to the next-to-lightest SUSY particle (NLSP). We consider the 
case when the lightest neutralino (x?) is the NLSP [6 1, and de- 
cays promptly with the dominant branching fraction yielding 
a photon and an essentially massless gravitino (%f — > Gy) |7|. 
The two gravitinos escape detection, resulting in the final state 
YY+ $t +X, where X denotes leptons and jets produced in the 
decay chains |[8l. 

In UED models, extra spatial dimensions are predicted that 
are accessible to all SM fields. We consider the case of a single 
UED that is compactified with radius R c , resulting in a tower 
of states for each SM field, called Kaluza-Klein (KK) exci- 
tations, with the masses of these states separated by R~ 1 . At 
the Tevatron, the UED process with the largest cross section 
would be the production of pairs of first level KK quarks J5), 
followed by decay chains to the lightest KK particle (LKP), 
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the KK photon (y*). If additional larger extra dimensions also 
exist that are only accessible to gravity, the LKP is able to 
decay promptly through gravitational interactions to a photon 
and a graviton (y* — » Gy) ifTOlfTTI . The two gravitons escape 
detection, resulting in the final state yy+ fx + X. 

Searches for BSM physics in yy+ fx +X events have been 
performed at the CERN e+e~ Collider (LEP) JT2), and at 
the Tevatron in Run I [13] and Run II lfl4Ul7l . This anal- 
ysis uses similar methods to those adopted in Ref. 1171 . a 
six times larger dataset, and improved photon identification 
criteria utilizing a neural network (NN) discriminant recently 
employed in other analyses fl8l . The larger dataset has sub- 
stantially increased the search sensitivity, and has allowed 
an improved formulation of the data-derived SM background 
prediction. The background prediction, including the assess- 
ment of systematic uncertainties, was developed using only 
the fx < 50 GeV region of the yy sample. Once finalized, 
the events with fx > 50 GeV were included in evaluating the 
consistency with the SM prediction and the sensitivity to the 
signal models. In addition to substantially improved limits on 
the GMSB model, this Letter also presents the first limits on 
the UED model with gravitational decays. 

The DO detector Hl9l consists of an inner tracker, a liquid- 
argon/uranium calorimeter, and a muon spectrometer. The 
tracking system is comprised of a silicon microstrip tracker 
(SMT) and a central fiber tracker (CFT), both located within a 
2 T superconducting solenoidal magnet. A central calorime- 
ter (CC) covers pseudorapidities |t|| < 1.1, and two endcap 
calorimeters (EC) extend the coverage to |t|| < 4.2, where 
T| = — ln[tan(0/2)], and 9 is the polar angle with respect to the 
proton beam direction. The electromagnetic (EM) section of 
the calorimeter is segmented in four longitudinal layers (EM/, 
i= 1,4) with transverse segmentation Ar| x A(|) = 0.1 x 0.1 (<]) 
is the azimuthal angle), except in EM3 where it is 0.05 x 0.05. 
A central preshower detector (CPS) utilizing several layers of 
scintillating strips, positioned between the solenoid coil and 
CC, provides a precise measurement of EM shower position. 
The trajectory of photon candidates is reconstructed by com- 
bining the four EM-layer and CPS measurements ifTTl . 

The data analyzed were collected with single EM trig- 
gers and correspond to an integrated luminosity of 6.3 ± 
0.4 fbT 1 ll20l . Events containing identified calorimeter noise 
patterns which could bias the fx distribution are removed. 
Diphoton candidate events are selected by requiring at least 
two photon candidates with transverse energy Ej > 25 GeV 
identified in the CC. Photon candidates are selected from 
EM clusters reconstructed within a cone of radius %_ = 
v/(Ar|) 2 + (A(j)) 2 = 0.2 by requiring (i) > 95% of the clus- 
ter energy be deposited in the EM layers, (ii) the calorime- 
ter isolation variable / = [£■,„, (0.4) — Eem(0.2)] /Eem(0-2) be 
less than 0.10, where E tot (2Q (E EM (2Q) is the total (EM) en- 
ergy in a cone of radius , (iii) the shower width in EM3 
be consistent with an EM shower, (iv) the scalar sum of the 
transverse momentum (px) of tracks originating from the pp 
collision vertex (PV) in a 0.05 < %_< 0.4 annulus about the 
cluster centroid be less than 2 GeV, and (v) the cluster not 



be spatially matched to a reconstructed track or a significant 
density of SMT and CFT hits IfTTl . Further rejection of jets 
misidentified as photons is achieved with a requirement on 
the NN discriminant, trained using a set of track, CPS, and 
calorimeter based variables lfl8l . 

Electrons satisfy the same requirements as photons, with 
the exception of the track veto (item v). Jets are recon- 
structed with the iterative midpoint algorithm [21] with cone 
size %_ — 0.5. The fx is determined using calorimeter energy 
depositions with |t|| < 4. Corrections are applied to fx to ca l- 
ibrate energy from EM objects and jets, and to account for the 
px of muons. There are on average several pp interactions per 
crossing of the beams. The correct PV is identified in ss 98% 
of signal events for the benchmark models. The photon tra- 
jectories must indicate that the candidates originate at the PV. 
This requirement is to ensure an accurate calculation of trans- 
verse energy in background events in which the correct PV 
is less efficiently identified, to suppress non-collision events, 
and measured to be s» 86% efficient using a Z(— > ee,p.p) +y 
data sample. To reduce the number of events with signifi- 
cantly mismeasured fx > events are rejected if the difference 
in azimuthal angle (A(j)) between the highest Ex jet (if present) 
and fx is greater than 2.5 radians, or if A(j) between either pho- 
ton and fx is less than 0.2 radians. A total of 7934 yy candi- 
date events satisfy these criteria. 

SM background events in the yy sample are categorized as 
arising from instrumental fx sources (SM yy, y+jet, multi- 
jet) and genuine fx sources (Wj, W+jet, W/Z + yy). All 
backgrounds are measured using data control samples, with 
the exception of small contributions from W/Z + yy events, 
which are estimated using Monte Carlo (MC) simulation. 

Instrumental fx is a result of energy mismeasurement in an 
otherwise Ex balanced event. Instrumental fx sources in the 
yy sample are separated into contributions from SM yy events, 
and events with at least one photon candidate originating from 
a misidentified jet (misID-jet), i.e., y+jet and multi-jet events. 
The difference in energy resolution for real photons and fakes 
from misidentified jets results in a difference in the shape of 
the fx distribution between the two categories. 

The fx shape in SM yy events is modeled using a dielec- 
tron (ee) data sample predominantly composed of Z — » ee 
events. The ee sample satisfies the same kinematic require- 
ments as the yy sample, with the exception that the ee invari- 
ant mass is restricted to an interval about the Z boson peak 
to reduce genuine fx contributions (e.g., W+jet, di-boson, 
and it events). The fx distribution in ee events is compared 
with shapes in Z — > ee and SM yy MC events generated with 
PYTHIA J22). These MC samples, and all others used in this 
Letter, were processed with full GEANT l23l detector simula- 
tion and standard reconstruction algorithms. Kinematic dif- 
ferences between the Z — > ee and SM yy processes are verified 
with MC to have a negligible impact on the fx shape. The 
Z — > ee MC accurately models ee data for fx values below 
fx ~ 35 GeV. Above this value, a more pronounced tail is 
observed in ee data. The tail in data reflects both mismeasure- 
ments not modeled in MC, and a small residual presence of 
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genuine events in the ee sample. The average of the data 
and MC shapes is used to model the in SM yy events for 
values of > 35 GeV, and the data-only and MC-only ex- 
tremes are used to define a systematic uncertainty on the this 
shape. 

The fir shape in misID-jet events is modeled with a data 
sample satisfying the same requirements as the yy sample with 
the exception that at least one of the photon candidates fails 
the NN requirement. Additionally, photon identification re- 
quirements (iii) and (iv) are loosened to reduce the statistical 
uncertainty on the fj shape. A systematic uncertainty on the 
fj shape in events with misidentified jets is obtained by vary- 
ing the photon identification criteria. 

The instrumental fj background estimate is normalized 
such that the number of events with < 10 GeV is equal 
to that in the yy sample. The relative contribution of SM yy 
and misID-jet background events is determined by a fit to the 
yy sample f-p distribution for fj < 20 GeV. The fit accounts 
for the small contribution of SM background with genuine 
fx m the fit region, and is verified to be insensitive to signal 
contributions for benchmark model cross sections relevant to 
this analysis. The SM yy contribution to the yy sample over 
the full grange is (41 ± 17)%. A systematic uncertainty ac- 
counts for changes in the shape of the predicted instrumental 
fj distribution arising from the uncertainty in the determina- 
tion of the SM yy contribution. 

SM background with genuine f-p arises from real SM yy+ 
$t +X events and from events with an electron misidenti- 
fied as a photon (misID-ele). The misID-ele contribution is 
derived using an ey data sample, composed primarily of in- 
strumental fj sources for fir < 20 GeV and W(— > ev)y and 
W(— > ev) + jet events at higher Rvalues. The instrumental 
fx sources are modeled with the previously introduced ee and 
misID-jet shapes, respectively. The Z->ee normalization 
is determined by fitting the Z boson peak in the ey invariant 
mass distribution, and the multi-jet fir shape is normalized to 
provide the remaining contribution in the fj- < 10 GeV re- 
gion. The presence of real fj contributions in the ey sample 
is seen as an excess of events with high values above the 
predicted contributions from instrumental sources. This ex- 
cess is well described by Wy and W + jet events. The expected 
W boson peak is observed in the transverse mass distribution 
of ey sample events with fj > 30 GeV. The normalization 
of the W + jet contribution is determined from a comparison 
of the data photon NN shape with MC real and fake photon 
NN shapes [18] in this $y region. The remaining contribu- 
tion is in good agreement with PYTHIA Wy production after 
applying a next-to-leading order (NLO) QCD correction 1241 
and an additional +15% scaling factor accounting for QED fi- 
nal state radiation (FSR) in inclusive W production. The FSR 
component l25ll is determined with data using the A^.(e,y) 
distribution. The predicted misID-ele contribution to the yy 
sample equals the excess of high fj- events in the ey sample, 
scaled by / e _>.y/(l — /<>-w), where / e _yy — 0.020 ± 0.005 de- 
notes the rate at which an electron fakes a photon satisfying 
the selection criteria, as measured with Z — > ee data. 
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FIG. 1 : f.j distribution in the yy sample shown with statistical uncer- 
tainty and expected SM background from events with a misidentified 
jet, a misidentified electron, W/Z + yy events, and SM yy events. The 
expected $j distribution in the presence of GMSB and UED events 
is also displayed for example values of A and R~ 1 , respectively. 



Real SM diphoton events with large genuine originate 
from W/Z + yy processes. This small background contribu- 
tion is estimated with MC using MADGRAPH (26). Events 
with inclusive W and Z boson decay modes are simulated, 
with W — > Iv (I = e,fj,z) and Z — > vv providing the largest 
genuine contribution. A total of 1.6 ±0.1 W+yy events 
and 3.8 ± 0.3 Z + yy events are estimated to be present in the 
yy sample. Figure [T] displays the yy sample fj distribution, 
which is in good agreement with the SM prediction over the 
full ^7- range. Table [I] provides the observed number of yy 
sample events and the SM prediction in three regions. 

We determine the sensitivity to the GMSB scenario using 
a set of values, termed SPS8 11271 . for the model parameters. 
In this set the scale A is unconstrained, M mes — 2A, N mes = 1, 
tanp = 15, and /j > 1271 . The masses and decay widths 
of SUSY particles are calculated with susyhit 1.3 J28) and 
used to generate PYTHIA MC events. The event selection ef- 
ficiency is 0.17 ± 0.02 at A = 120 TeV, and does not differ 
significantly for other A values studied. The NLO production 
cross section is calculated with PROSPINO 2.1 Q. The ex- 
pected $ T distribution for the SM and GMSB at A = 120 TeV 
is depicted in Figure [T] The number of expected GMSB 
events in three fir regions is listed in Table|l]for A = 100 and 
120 TeV. 

We consider the UED model as implemented in PYTHIA 
6.421 [29 1, leaving R~ l unconstrained and setting AR C = 
20, where A is the cutoff scale for radiative corrections to 
KK masses. This UED model is implemented in a higher 
(4+Af) dimensional space, where R~ 1 is much larger than 
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TABLE I: Observed number of yy sample events, predicted background from instrumental fj sources (SM yy, Y+jet, QCD multi-jet) and 
genuine IjLj sources (Wy. W +jet, W/Z + yy), and total predicted SM background, in three fj intervals. The expected number of GMSB and 
UED signal events is listed for two A and R~ values, respectively. The total uncertainty on the SM background and expected signal is given. 



fx Interval, 


Observed 


SM Background Events 






Expected Signal Events 




GeV 


Events 


Instr. fir 


Genuine $j 


Total 


GMSB 


GMSB UED 


UED 












A = 100 TeV 


A = 120 TeV R- 1 = 420 GeV 


R- 1 =460 GeV 


35-50 


18 


9.6±1.9 


2.3±0.5 


11.9±2.0 


1.8±0.1 


0.3±0.1 1.4±0.1 


0.3±0.1 


50-75 


3 


3.5±0.8 


1.5±0.3 


5.0±0.9 


4.1 ±0.3 


0.8±0.1 2.9±0.2 


0.6±0.1 


>75 


1 


1.1 ±0.4 


0.8±0.1 


1.9±0.4 


14.3 ±1.1 


4.4 ±0.4 24.7 ±2.0 


6.4±0.5 



that of the N compact extra dimensions accessible to grav- 
ity, inducing KK particle decays through gravitational inter- 
actions. We choose N = 6 and a fundamental Planck scale 
Mo = 5 TeV, such that only the y* — > Gy decay occurs with 
appreciable branching fraction ifTTll . The event selection ef- 
ficiency is 0.19 ±0.02 at R~ l = 460 GeV, and does not dif- 
fer significantly for other R~ 1 values studied. The expected 
$ T distribution for the SM and UED at R~ l = 460 TeV is 
depicted in Figure [T] The number of expected UED events 
in three regions is listed in Table [i] for R~ 1 = 420 and 
460 GeV. 

Systematic uncertainties for sources of instrumental ft? are 
attributed to the uncertainty of the fj shape in SM yy and 
misID-jet events, and their relative normalization. An un- 
certainty in the shape of the distribution for the misID- 
ele contribution arises from the uncertainty in the Z — » ee 
contribution to the ey sample, and a 25% misID-ele normal- 
ization uncertainty results from the / e -yy uncertainty. Sys- 
tematic uncertainties in the contributions estimated with MC 
arise from the integrated luminosity (6.1%), trigger efficiency 
(2%), and photon identification (3% per photon) and trajecto- 
ries (3%) efficiencies. Uncertainty in parton distribution func- 
tions (PDF) [30 1 yield systematic uncertainties of up to 5% 
and 20% in the production rate of GMSB and UED events, 
respectively. 

No evidence for BSM physics is observed in the yy sample 
fj distribution and limits on the benchmark models are de- 
rived using a Poisson log-likelihood ratio test [31] incorporat- 
ing the full fj- distribution. Pseudo-experiments are generated 
according to the background-only and signal plus background 
hypotheses, and account for statistical uncertainty on the ex- 
pected number of events and systematic uncertainties. The 
cross section limit is evaluated using the CL V modified fre- 
quentist approach ||3TI . Figure [2] shows the predicted GMSB 
and UED cross section with PDF uncertainty, and 95% C.L. 
cross section exclusion limit, as functions of A and R~ l , re- 
spectively. For GMSB, the NLO cross section uncertainty is 
small compared to the PDF uncertainty. The UED NLO cross 
section has not yet been computed. 

In conclusion, we have presented a search for physics be- 
yond the standard model in the yy + fx +X final state at 
the Tevatron. The observed fj distribution is consistent with 
the SM expectation and limits on two benchmark models are 
derived. In the SPS8 GMSB model, values of the effective 



SUSY breaking scale A < 124 TeV are excluded at 95% C.L. 

The limit excludes m o < 175 GeV, representing improve- 

^i 

ments of 50 GeV [ 17 1 and 26 GeV [ 15 1 with respect to previ- 
ous measurements. Additionally, the first assessment is made 
of the sensitivity to the UED model with KK particle decays 
induced by gravitational interactions, excluding values of the 
compactification radius R~ l < All GeV at 95% C.L. 
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